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A recently derived method [R. D. Rohrmann and A. Santos, Phys. Rev. E. 76, 051202 (2007)] 

to obtain the exact solution of the Percus-Yevick equation for a fluid of hard spheres in (odd) d 

dimensions is used to investigate the convergence properties of the resulting virial series. This is 

Jh . done both for the virial and compressibility routes, in which the virial coefficients Bj are expressed 

^T^' in terms of the solution of a set of (d — l)/2 coupled algebraic equations which become nonlinear 

'^^ , for d > 5. Results have been derived up to d = 13. A confirmation of the alternating character of 

f0y. • the series for d > 5, due to the existence of a branch point on the negative real axis, is found and 

-isj ' the radius of convergence is explicitly determined for each dimension. The resulting scaled density 

per dimension 2ri^''', where rj is the packing fraction, is wholly consistent with the limiting value 

of 1 for d —> oo. Finally, the values for Bj predicted by the virial and compressibility routes in 

f"1~! ' the Percus-Yevick approximation are compared with the known exact values [N. Clisby and B. M. 

O ! McCoy, J. Stat. Phys. 122, 15 (2006)]. 

C/3 ' 

2 I I- INTRODUCTION equation of state both from the virial and compressibility 

routes in a rather straightforward procedure. Also very 
^ Interest in studying fluids of hard spheres recently, Adda-Bedia et al^ have been able to reduce 

S in d dimensions goes back at least four the PY equation for hard disks (d = 2) to a set of sim- 
O ' j^^. 3^ ji; ji.2.3.4.5.6.7.8.9.io. 11. 12. 13. 14. 15. 16. 17.18.19.20.21. 22.23.24.25.26.rilea ntegral equations which they then solve numerically. 
O: and has recently experienced a new An interesting aspect of this work is their claim that the 

Ijpp^ 29.30.31.32.33.34.35.36.37.38.39.40.41.42.43.44.45.46.47.48.49.50.51.3aieSlia(i5i;Bgty?'l3e.$£eKe ralized tO any evCn dimension. 

The evaluation of virial coefficients, the derivation of the Due to the limited availability of virial coefficients and 

equation of state for these systems, or the determination to the fact that the virial series for high densities relevant 

of the radius of convergence of the virial series arc to the fluid phase seems to be in general either divergent 

J^ ■ among the issues that have been examined, but many or slowly convergent, various scries convergence accclcr- 

questions related to these issues and others are still ating methods, such as Fade or Levin approximants, have 

open. Providing answers to these questions may shed been used to derive approximate equations of state for d- 

^ ', light on similar issues related to real fluids and therefore dimensional hard-sphere fluids. However, the radius of 

^ ■ efforts in this direction are called for. convergence of the virial series is not known in general, 

O ; The Percus-Yevick (PY) theory^ is one of the classic and hence the usefulness of such approximate equations 

L; ■ integral-equation approximations of liquid-state theory,^ of state is limited by the uncertainty of its range of appli- 

.^ ; Apart from yielding the correct expression for the radial cability For d < 3 all known virial coefficients turn out 

V '. distribution function (rdf ) to first order in the density to be positive but since negative virial coefficients appear 

CI ■ (and hence also the correct second and third virial co- ^^r d > 4, the question arises as to whether some higher 

. F^. ! efficients), its key role in the theory of simple liquids order virial coefficients both for d = 2 and d = 3 might 

was recognized very early because the resulting integral eventually become negative. 

equation is analytically solvable for the important case of The main purpose of this paper is to examine the ques- 

the hard-sphere fluid {d = 3)i^i^ Further, the approx- tion of convergence of the virial series for fluids of d- 

imation provides the exact rdf (although not the exact dimensional hard spheres by looking at approximate thc- 

cavity function inside the core^) for hard rods (d = 1) orics where both the virial coefficients and the equation of 

and has been proven to yield exactly solvable equations state are known. We will use the procedure introduced in 

in odd d dimensions j^>^ In fact, explicit analytical solu- Ref. l58l to derive the equation of state of the system and 

tions for d = 5 and d = 7 have been obtaine d^'^'^^i^^ the values for the virial coefficients taking both the virial 

and, rather recently, numerical solutions for d = 9 and and compressibility routes, all within the PY approxima- 

d = 11 have been reportedi^ These latter have been tion. In particular, we will investigate the convergence 

derived using an alternative method to the one origi- properties of the virial series of the equations of state 

nally employed by Leutheusser— which, amongst other stemming out of both routes for d ~ 5,7,9,11, and 13 

things, allows one to obtain the virial coefficients and the and will show that the radius of convergence is related 



to a branch point on the negative real axis. Moreover, 
we will compare the PY virial coefficients with the ex- 
act values available in the literature. This comparison 
suggests that, as the dimensionality increases, the true 
radius of convergence tends to the value predicted by the 
PY theory. We will also examine the performance of the 
compressibility factors obtained from both routes with 
the corresponding simulation data and will show that the 
virial series truncated just before the first negative term 
provides an excellent approximation to the equation of 
state of the fluid phase. 

The paper is organized as follows. In the next section 
and in order to make the development self-contained, we 
provide a brief description of the so-called Rational Func- 
tion Approximation (RFA) approach leading to the PY 
approximation for fluids of hard hypersphcres in odd di- 
mensions. Section IIIII presents the numerical results of 
our calculations (together with a comparison with sim- 
ulation data) for the PY compressibility factors derived 
via the virial and compressibility routes, respectively, as 
well as the analysis of the behavior of the PY virial coef- 
ficients obtained from the same routes to the equation of 
state both with respect to the convergence properties of 
the virial series and with respect to the exact known val- 
ues. The paper is closed in Sect. |IV] with some discussion 
and concluding remarks. 



II. THE PY THEORY 

In this section we provide an outline of the main steps 
leading to the PY approximation for the thermodynamic 
and structural properties of hard hypersphere fluids in 
odd dimensions. Instead of following the original deriva- 
tion by Leutheusser^ we will use the RFA method intro- 
duced in Ref. [58i which the reader is urged to consult for 
further details. 



A. General background 

We begin by recalling that the static structure factor 
of a d-dimensional hard-sphere fluid can be written as 






G{-tk)] , (2.1) 



where p is the particle number density and G{s) is a 
Laplace-space functional given by^ 
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Here g{r) is the rdf and 9n{t) is the reverse Bessel polyno- 
mial of degree n = (d — 3)/2. From these structural func- 
tions one may obtain thermodynamic quantities such as 
the compressibility factor Z = p/ pksT and the isother- 
mal susceptibility x = kBT{dp/dp)T; both defined in 



terms of the temperature T, the pressure p and the num- 
ber density p {ks being the Boltzmann constant). The 
actual relationships read 



and 



1 + 2^-i775(cT+) 



X = ^(0), 



(2.3) 
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where a is the diameter of the particles and 77 is the 
packing fraction given by 
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ri = Vdpa 



Vd = 



(7r/4)'^/^ 
r(l + d/2)' 



(2.5) 



Vd being the volume of a d-sphere of unit diameter. 
Henceforth without loss of generality we will set a = 1. 

In hard-particle systems the temperature does not play 
any relevant role on the structural properties here intro- 
duced. Moreover, the thermodynamic state of such ffiiids 
can be characterized by a variable alone, say the density 
or the packing fraction. Therefore, taking into account 
the thermodynamic relation 
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the so-called virial equation of state (|2.3p and the com- 
pressibility equation of state (|2.4p provide two alterna- 
tive routes for obtaining the compressibility factor Z(r]) 
of a hard d-sphere fluid. Since all well-known theoretical 
methods to obtain structural functions give approximate 
results (with the exception of the exact solution for the 
one-dimensional case d = 1), one typically obtains two 
approximate solutions, the compressibility factor from 
the virial route 



Z„(r?) = l + 2'^-i77g,(,7), 



(2.7) 



with (?c(^) = 5(1^)7 Sijid the compressibility factor via 
the compressibility route 



^c(7?) = / dxx ^irix), 
Jo 



(2.8) 



with xiv) given by Eq. (|2.4p . 

For hard-hypersphere fluids in arbitrary odd d dimen- 
sions, the RFA approach provides an analytical repre- 
sentation of the function G{s) defined by Eq. (|2.2p and 
related to the structure factor (|2.ip . In the simplest im- 
plementation of the RFA approach, which is the so-called 
standard RFA and coincides with the PY approximation, 
the function G{s) can be written in the explicit form^ 
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s 1 + A?;^" aj^d-i(s) 



A=(-l)('^-i'/V-id!!, 
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(2.10) 
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The coefficients a-,- are in general functions of the density. 
Specifically, ao = (d — 2)!! and the quantities aj with 
j = 1, . . . , n + 1 are solutions of the following closed set 
of n + 1 equations: 
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E 72j£'2(m-j)-i =0, < m < n, 



j=o 



with the boundary condition 

{2n + 2-jy. 



i,,=o 2"+i-J(n + l-j)!j!' 



(2.12) 



(2.13) 



The coefficients Di are linear combinations of the {aj} 
given by Dq = 1 and 
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I > 1, 



and the quantities 72m with < m < n are given in terms 
of the coefficients aj by means of the recursion relation 
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E 72j-D2(m-j) , < m < n. 



j=o 



(2.15) 
Here we have adopted the conventions Uj = if j > ?i + 1 
and J2T=o • • • = if TO < 0. In summary, when the {72m} 
obtained from Eq. (|2.15p are inserted into Eq. (|2.12p . 
ans use is made of Eq. (|2.14p . one gets a closed set of 
n + 1 = (d — l)/2 algebraic equations for ai, 02, . . . , 
a„-|-i. The number of mathematical solutions (including 
complex ones) is 2" = 2^'^"'^^/^ and the physical solution 
is obtained as the one yielding the correct low density 
limit given by Eq. ((2T3)) . 

The contact value of the rdf and the isothermal sus- 
ceptibility obtained by the PY theory (or, equivalently, 
by the standard RFA method) are 



.9c (j?) = On+l 



"+1 ,' 1 \d-j 
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B. Virial expansions 

The virial expansion of the equation of state is an ex- 
pansion in powers of the density or the packing fraction, 



z{v) 



J=2 



nJ-1 
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and its range of validity is limited by the convergence 
properties of the series. Notice that 
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Since gdi] = 0) = 1, Eq. dm yields 



(2.19) 



(2.20) 



A more elaborated analysis of the virial equation shows 
that 63 = 2'^^^^yi{l)/yi{0), where yi{r) represents the 
overlap volume of two hyperspheres of unit radius and 
centers separated by a distance r. In particular, for d = 
odd, one finds^ 
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(2.21) 
The application of the PY theory yields two virial expan- 
sions, one for Z^ir]) and another for Zc(r\). The virial 
coefficients in the virial route are given by 



6(-) = 
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(j-2)! ** ■^"^"i';=o^ 
whereas in the compressibility route one has 
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with gcW) ^-iid xW> given by Eqs. (|2.16[) and (|2.17p . 
respectively, and where 9^ denotes the derivative with 
respect to 77. In practice, what one does is to solve Eq. 
(|2.12p in a recursive way for the coefficients in the den- 
sity expansion of oi, . . . , a„-|_i. The solutions are exact 
rational numbers and from them one gets \y" and h!y 
also as exact rational numbers. 



Cif]) = 1 - 2\{d - 2y\i^[Dd + Y. l2jDd-2-2,) , (2.17) 



j=0 



with the factors Di and 72m given by Eqs. (|2.14p and 
(|2.15[) . respectively. With these results introduced into 
Eqs. (|2.7p and (|2.8p one obtains the compressibility fac- 
tors Zy(r]) and Zc{r]) as derived within the PY theory 
from the virial and compressibility routes, respectively. 



III. RESULTS 

In this section we present the results that follow from 
the previous derivations. Three aspects will be analyzed. 
We first deal with the virial coefficients. Then, we exam- 
ine the issue of the convergence properties of the virial 
series and finally we compare the resulting compressibil- 
ity factors with simulation data. 



TABLE I: Exact and PY values for the virial coefficients 
fej/b^~^ for 4 < J < 10 and several dimensionalities. The 



ues for hard disks (d = 


- 2) are obtained from Ref. [63- 


Exact 


PY (v) PY (c) 


d bi/hl 


CV&i ^Vl^l 



2 
3 
5 
7 
9 
11 
13 
15 
d 
2 
3 
5 
7 
9 
11 
13 
15 
d 
2 
3 
5 
7 
9 
11 
13 
15 
d 
2 
3 
5 
7 
9 



0.53223180 
0.2869495 
0.07597248 
0.00986495 
-0.008581 
-0.011334 
-0.009523 
-0.006934 



bs/fet 



0.33355604 
0.110252 
0.0129551 
0.0070724 
0.007439 
0.006176 
0.004307 
0.002705 



bs/bi 



0.1988425 
0.03888198 
0.0009815 
-0.0035121 
-0.004794 
-0.00395 
-0.002580 
-0.001472 



h^jhl 



0.1148728 
0.01302354 
0.0004162 
0.0025386 
0.003716 



bs/b^ 



0.0649930 

0.0041832 

-0.0001120 

-0.0019937 

-0.003222 



bg/bj 



0.0362193 
0.0013094 
0.0000747 
0.0016869 
0.003029 



bio/bj 



0.0199537 
0.0004035 
-0.0000492 
-0.001514 
-0.00306 



0.5008 
0.25 
0.04785156 
-0.007499695 
-0.018590778 
-0.016933400 
-0.012604546 
-0.008614616 



h\rih% 



0.2948 
0.0859375 
0.01294708 
0.01235023 
0.01172436 
0.00873948 
0.00566100 
0.00337395 



CV&i 



0.1667 
0.02734375 
-0.00286007 
-0.00817700 
-0.00820391 
-0.00583856 
-0.00347797 
-0.00186395 



hVih\ 



0.09194 
0.00830078 
0.00245853 
0.00655314 
0.00676771 



hVlhl 



0.04988 

0.00244140 

-0.00173717 

-0.00576280 

-0.00621239 



h\rih\ 



0.02677 
0.00070190 
0.00134572 
0.00542528 
0.00615464 






0.01424 

0.00019836 

-0.00108546 

-0.00537523 

-0.00645733 



0.5378 

0.296875 

0.08905029 

0.021550816 

0.000373515 

-0.005004249 

-0.005256304 

-0.004146893 



bTlh\ 



0.3433 
0.12109375 
0.01638031 
0.00511681 
0.00404299 
0.00320718 
0.00221818 
0.00138891 



r/bi 



0.2090 
0.04492188 
0.00145829 
-0.00186533 
-0.00229461 
-0.00179806 
-0.00113940 
-0.00063735 



T^ 



/hi 



0.1233 

0.015625 

0.00063127 

0.00138425 

0.00164806 



bt'/hl 



0.07107 

0.00518799 

-0.00022064 

-0.00107878 

-0.00133333 



1^ 



/hi 



0.04025 
0.00166321 
0.00017676 
0.00091218 
0.00117921 



!,( = ) Ili) 
Old I "2 



0.02248 

0.00051880 

-0.00012788 

-0.00081922 

-0.00111652 



Virial coefficients 



Because the PY theory is exact to first order in density, 
the virial coefficients 62 and 63 [Eqs. ^^TMf and ^I^] 
are exactly reproduced by both routes. Higher virial co- 
efficients are, however, different for each route, as shown 
in Table H] for 4 < j < 10, where the known exact values 
arc also included for the sake of comparison. The normal- 
ized differences {h/" —bj)/\bj\ between the approximate 

virial coefficients \n and of and the exact values^i^i^ 
bj are shown in Fig. [1] as functions of the space dimension 
for 4 < j < 10. As one can see from the results in this 
figure, all known (exact) coefficients in d = 2 and rf = 3 



dc) 



lie between b\ ' and b\ ' in the form b 



(v) 



M 



< h < bj 



{c) 



but 



J.J. . 3 

this is not so in higher dimensions. A transition behavior 

seems to take place at d = 5 since in that case one has 

bj < bf < 5^"^ for j = 7 and 9, while bj > bf^ > bf'' for 

j = 8 and 10. For d> 7, however, the trends seem to be 



6^"^ < b, < b'f'' for j 



(c) (v) 

even > 4 and by < bj < b- for 

7, 



J 



odd > 5. The top panels of Fig. [Ijseem to indicate 



that the relative deviations of the PY values with respect 
to the exact ones tend to decrease and stabilize with in- 

(c) 

creasing d, especially in the case of bj . Of course, a 
confirmation of all these trends would require the knowl- 
edge of the exact virial coefficients for higher orders j and 

(v) (c) 

higher dimensionalities d. The fact that 64 < 64 < 64 
for all d implies that Zy{r]) < Z[ri) < Zc{r]) for asymptot- 
ically low densities. On the other hand, since both Zy{r]) 
and Zc{ri) keep being finite for any 77 < 1^^ i.e., even 
for densities higher than the close-packing value Jjcp, it 
can be reasonably expected that Z[ri) > Zc{i]) > Zy{ij) 
beyond a certain density, although this possibly happens 
in the metastable region. A precursor of this effect might 



h(=) 



(v) 



be the relation b^' < 65 < 6^ -' for d>7. 

The assessment of the performance of the PY virial 
coefficients with respect to the exact results may also 
profit from a different representation of the data. This is 
shown in the four panels of Fig. [21 where we have plotted 
the ratio b2bj^i/bj as a function of j for 3 < j < 20 



dc) 



(with both bj ' and b- ' and with the values from Refs. 
llmiii) for d = 2, 3, 5, 7, and 9. As we will discuss 
below, the magnitude of this ratio is related to the radius 
of convergence of the virial series. In this instance, the 
regularity of the results for the lower j's observed for 
d = 2 and d = 3 is lost in higher dimensions. It is 
interesting to note that for d = 7 and 9 the exact values 
of b2bj-i/bj for the higher j's (7 < j < 10) lie very 
close to the PY values. Whether this good agreement is 
accidental or not cannot be assessed before exact values 
of bj for j > 11 and/or d > 11 are known. 



B. Radius of convergence of the virial series 



Now we turn to the question of the convergence prop- 
erties of the virial series. The radius of convergence of the 
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FIG. 1: (Color online) Normalized differences between the 
virial coefficients 6'" (triangles) and b '^ (circles) and the ex- 
act value a^^'^^'^^ bj, as functions of the space dimension for 
4 < j < 10. The data for d — 2 have been obtained from Ref. 
mi . The lines have been drawn to guide the eye. 
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FIG. 2: (Color online) Ratio b2fej-i/6j as a function of j 
for 3 < i < 20 and d = 2, 3, 5, 7, and 9; fe*"' (triangles), 

by (circles), and exact value^'^'^ bj (squares). The dash- 
dotted line is the limiting value of this ratio for both PY virial 
coefficients when j — > oo. The insets in the cases d = 5, 7 and 
9 show in more detail how the limiting value is approached. 



TABLE II: The radius of convergence 770 of the virial series in 
the PY theory, and the associated quantities &2?7o and po = 



2r;Q' for fluids at d = 


= 3, 5, 7, 9,11, and 13. 




d rjo 


biVo 


PO 



3 
5 
7 
9 
11 
13 



0.056624327 
0.010886779 
0.002394546 
0.000557016 
0.000133287 



0.90598923 

0.69675386 

0.6130038 

0.570384 

0.545943 



1.58740105 
1.12623911 
1.04854495 
1.02289857 
1.01204491 
1.00678503 



virial scries for each dimension d, r/o = liiiij^oc \bj-i/bj\, 
is determined by the modulus of the singularity of Z(ri) 
closest to the origin in the complex rj plane. In order 
to inhibit the influence of d on the characteristic density 
values, we will sometimes choose 6277 (rather than 77 or 
pa"^) to measure the density. In Fig. [3] we display the 
PY compressibility factor Zyirf) as a function of &2'7 for 
d = 5, 7, 9, and 11. In the figure, apart from showing 
Zy{ri) in the physical domain of positive densities (thick 
solid line, shaded region), we have provided its analytic 
continuation to negative values of 77. The thin solid line 
shows such a continuation. It turns out that there ex- 
ists a certain negative value rj = —t]q such that Z„{vi) 
keeps being real in the interval —770 < 77 < 0. How- 
ever, at 77 = —r/o Zy{r]) merges with an unphysical root 
(dashed line) and both roots become a pair of complex 
conjugates for 77 < —770 (the two dotted lines represent 
their imaginary parts). This shows that Z^, (77), as well 
as Zc{ri), possesses a branch point at 77 = —770. This is 
the singularity on the real axis closest to the origin. If 
no other singularity lying in the complex plane is closer 
to the origin, then 770 is the radius of convergence of the 
series. Figure |4] provides the radius of convergence of the 
virial series for Zv{t]) and Zc(ji), this time by represent- 
ing again b2bj'Zi /uf as a function of j for 3 < j < 150 
and fi = 3, 5, 7, 9, 11, and 13. As is well-known, the ra- 
dius of convergence predicted by the PY theory is 770 = 1 
for d = 3. On the other hand, for d > 5 the radius of 
convergence is 770 < 1 and coincides with the value 770 
corresponding to the branch point on the negative real 
axis shown in Fig. [31 The values of 770, 627701 and of the 
scaled density per dimension^ po = 277q are shown in 
Table HIl for odd dimensions in the interval 3 < d < 13. 
It can be observed that the values of po are consistent 
with the limit po — ^ 1 as d — > 00 conjectured by Frisch 
and Percusii^ This agreement, along with the behavior 
observed in Fig. ^ supports the reliability of the radius 
of convergence predicted by the PY theory, at least for 
high dimensions. 



C. Compressibility factors 

We now turn to the compressibility factor. In Fig. 
[5] we present a comparison between the PY values for 
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FIG. 3: (Color online) Compressibility factor Z^(r\) in the 
physical domain of positive densities (thick solid line and 
shaded region) and its analytic continuation for 77 < (thin 
solid line) as a function of 6277 for d = 5, 7, 9, and 11. At 
the negative value 77 = —770 Zv{rj) merges with an unphysi- 
cal real root (dashed line) and both roots become a pair of 
complex conjugates for 77 < —770 whose imaginary parts are 
represented by the two dotted lines. The boxed numbers are 
the corresponding values of —62770. 



Zyirf) and Zc^r]) and some of the available simulation 
datai^i2id9i5ii53j54^^^£7 foj. yarious dimensions. It fol- 
lows from this figure that, as d increases, Zc{r]) becomes 
a rather accurate approximation for the simulation re- 
sults. Note, however, that the simulation data for d ~ 9 
are restricted to the density region where Zy{r]) ~ Zc{rj). 
Comparison between the density range in Fig. [5] and 
the values of r/o tabulated in Table [III shows that the good 
agreement between Zc and Z for d = 5 and d = 7 extends 
to 77 > 770, i.e., well beyond the radius of convergence 
of Zy and Zc- This might cast doubts on the practical 
usefulness of the virial coefficients to predict the equation 
of state of hard hypersphcres in the fluid region with 
77 > 770. However, as the following discussion shows, we 
have observed that this is not the case. Let us denote 
by j* the order of the virial coefficient preceding the first 
negative coefficient, i.e., bj > for j < j* but foj*+i < 0. 
Analogously one can define j* and j* in connection with 

&j" and &^^ , respectively. For instance, at d := 5 one has 
(see Table U j* = j* = 8 and j* = 6, while j* = j* = 6 
and j* = 4 at d = 7. Let us now define the truncated 
virial series 



J=2 



(3.1) 



FIG. 4: (Color online) Ratio fe2&j-i/&j as a function of j for 
3 < J < 150 and d = 3, 5, 7, 9, 11, and 13; 6'" (triangles) and 

b '^ (circles) . The dash-dotted line is the limiting value of this 
ratio for both PY virial coefficients when j — > 00 and from 
such value one can immediately get the radius of convergence 
of the virial series. 



with similar definitions of Z*{ij) and Z*{ij). Since the 
first neglected term is negative, it follows that, at least for 
low densities, Z*{r]), Z*{r]), and Z*{ri) are upper bounds 
oi Z{r]), Zy{rj), and Zc(rj), respectively. Furthermore, we 
have checked (for d = 5, 7, 9, 11, and 13) that Z*{r]) 
and Z*{r]) are excellent estimates of Zyirj) and Zc{ri), 
respectively, much better in the region r] > rjo than any 
other truncated virial series. Does this property extend 
to Z*{ri) as well? Comparison between Z*{r]) and sim- 
ulation data of Z{ri) for d = 5 and 7 indicate that this 
is indeed the case. In Fig. [6] we have represented the 
approximate compressibility factors Zy{r]), Z*{ri), and 
Z*{ri) for d = 5 and 7, and compared such approxima- 
tions with simulation values ^ii^i^ Again the superiority 
of the approximate PY compressibility factor obtained 
via the compressibility route over the one obtained via 
the virial route is apparent. More interesting, however, 
is the fact that if one considers the similar approxima- 
tion Z*(rf) (i.e., using the exact virial coefficients), the 
agreement between this approximation and the simula- 
tion data is strikingly good over the whole fluid phase 
region where these data are available, especially in the 
case d = 5. In fact, since as already pointed out, the 
PY theory leads to the exact (positive) 62 and 63 (ir- 
respective of whether one takes the virial or the com- 
pressibility routes), and since for d > 9 both the exact 

fourth virial coefficient and the PY 64 are negative, it 
turns out that for these dimensions j* = j* =3 and the 
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FIG. 5: (Color online) Compressibility factors Z^(r\) (dashed 
line) and Zc{rj) (continuous line) of the PY theory as functions 
of the packing fraction 77 for rf = 3, 5, 7, and 9, and simulation 



dataj 



21 .49.59.66.67 



(filled circles). 



FIG. 6: (Color online) Truncated approximations Z*{r]) 
(dashed line), Zc{r]) (continuous line), and Z* (rf) (dotted line) 
as functions of the packing fraction 77 for d = 5 and 7, and 



simulation dataf 



' (filled circles). 



truncated expansions Z*[ri) = Z*{ri) ^ 1 + 62'? + ^3*7^ 

(c) 

coincide. Furtermore, since the PY coefRcient b\ also 
becomes negative for d > 11, all three truncated expan- 
sions for the compressibility factor become the same for 
d> 11. 



IV. DISCUSSION 

The preceding results lend themselves to further 
consideration. One might have reasonably wondered 
whether the trend observed for all the known virial coef- 
ficients bj both in d = 2 and d = 3, namely the fact that 

they are bracketed by 6^" and 6^*^ , would remain valid for 
all virial coefficients in these dimensions and also hold for 
the higher dimensions. The following reasoning indicates 
that for any d > 1 there must be at least one virial coef- 
ficient that does not comply with the above trend. Since 
both Zy{ri) and Zc{vi) only diverge for rj = 1 and the true 
Z{ri) must have a divergence at the close-packing fraction 
77cp < 1, then at a smaller packing fraction (although pos- 
sibly in the mctastable ffuid region), Z{ri) must lie above 
Zyiji) and Zc(ri). This crossing should manifest itself in 
the behavior of the virial coefficients even if the radius of 
convergence of the virial series is smaller (or very much 
smaller) than the packing fraction at which the cross- 
ing occurs. According to our findings, the bracketing 

7, where one 



bj < bj < ¥p is lost already for 65 at d 



has b 



{c) 



< h < b, 



(v) 



, and it seems that this will happen 



for all higher odd j's but not for even ones. We also find 
that the relative deviations of the PY virial coefficients 
with respect to the true coefficients tend to stabilize with 
increasing dimensionality. 

Concerning the radius of convergence of the virial se- 
ries, two key aspects should be pointed out. The first one 
is that, as our present results clearly confirm, in the PY 
theory and for d > 5 the singularity closest to the origin is 
a branch point located on the negative real axis and cor- 
responds to a packing fraction 770 much smaller than the 
one delimiting the ffuid phase region. The second obser- 
vation is that the PY values for the ratio 626J-1/&J for the 
higher j's lie very close to the exact values. Further, the 
resulting scaled density per dimension po is wholly con- 
sistent with the limiting value of 1 for d —> 00. Hence, on 
the one hand the scenario of an alternating series for the 
true virial series by Clisby and McCoy^ (perhaps even 
for d = 3) is reinforced and, on the other, one may conjec- 
ture that its radius of convergence will also be close to the 
one in the PY theory. In connection with the latter point 
it is instructive to compare the PY values of the radius 
of convergence ryo with known bounds. In 1964, Lebowitz 
and Penrose^ derived the lower bound 62770 > 0.07238, 
which is clearly vcriffed by the values of Table |lll They 
also showed that no phase transitions exist for values of 
62?? less than 1/2(1 -I- e) ~ 0.135, which is again consis- 
tent with the PY predictions for 770 and the scenario^ 
of a leading singularity on the positive real axis at the 
freezing density rjf > tjq. Recently, Clisby and McCoy^ 
estimated the value 62^7 = 0.5 for the radius of conver- 



gence of the sum of the Ree-Hoovcr diagrams. It is quite 
interesting to notice that the values of 62770 predicted by 
the PY approximation (see Table |ll| seem to tend from 
above precisely to 0.5, i.e., lim^-^oo 62?7o = 0-5. 

As a final point, it is worth mentioning that the regu- 
larity observed for the ratio {Zy — Z)/{Zc — Zy) in d = 3 
and d = 5, which leads to the Carnahan-Starling equa- 
tion of stated in the former case and something similar 
in the latter (2^ seems not to hold for d ~ 7, even if the low 
density values (which are subjected to great uncertain- 
ties) are ignored. The recent simulation data of Bishop 
et al^ are also for low density and with error bars com- 
parable to Zc — Zy so that further conclusions on this 
matter are precluded at this stage. 
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